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Abstract 

We constructed an analog circuit generating fluctuations of which a proba- 
bility density function has power law tails. In the circuit fluctuations with 
an arbitrary exponent of the power law can be obtained by tuning the resis- 
tance. A theory of a differential equation with both multiplicative and addi- 
tive noises which describes the circuit is introduced. The circuit is composed 
of a noise generator, an analog multiplier and an integral circuit. Sequential 
outputs of the circuit are observed and their probability density function and 
autocorrelation coefficients are shown. It is found that correlation time of 

the autocorrelation coefficient is dependent on the power law exponent. 
* Electronic mail : aki@sawada.riec . tohoku .ac.jp 



1 Introduction 



Random multiplicative process (RMP) is attracting much interest recently 
as a new mechanism of generating power law fluctuations. The process is de- 
scribed by a Langevin system with both multiplicative and additive noises. It 
was proved that a PDF for its amplitude has power law tails when the noises 
are independent |l|, ^ ^. The continuous time version of RMP is described by 
a stochastic differential equation with both multiplicative and additive noises. 
It was shown that a stationary probability density function for its amplitude 
has power law tails when the noises are Gaussian white noise[|, ^, |^. This 
type of equations has been studied in chaos on-off intermittency0, motions 
of a polymer in turbulent fluid 0, nonlinear coupled oscillators ^ and price 
changes in economics []TU|. 



In this paper we propose an electrical generator of random fluctuations 
based on the theory of RMP for the power law fluctuations by introducing 
an analog electrical circuit consisted of operational ampliflers and an analog 
multiplier. 

The paper is organized as follows. In Sec. ^ we make a brief explanation 
of the theory of the continuous time RMP. In Sec. ^ we show the circuit 
diagram of the proposed generator having power law fluctuations and explain 
its mechanism. In Sec. ^ we show the data of sequential outputs, their 
probability distribution function (PDF) and their autocorrelation coefficients 
with discussion. Sec. |] is devoted for concluding remarks. 

2 Theory of RMP 

We here review a brief outline of the stochastic differential equation with 
both multiplicative and additive noises@, The stochastic differential 

equation is given by 



where v{t) is a dynamical variable, and z/(t) and ^(t) represent a multiplica- 
tive noise and an additive noise, respectively. Representing ensemble average 
by (...) we require = z?, and (^(t)) = 0. When z/(t) and ^(t) follow 

Gaussian white noises with their strengths characterized by Di, and D^, re- 
spectively, 

{[u{ti) - uMh) - u]) = 2DJ{ti-t2), (2) 
(e(ti)e(t2)) = 2D^S{h-t2), (3) 

the PDF for the dynamical variable v in eq. (|l]), p{v,t), is known to follow 
the generalized Fokker-Planck equation. 



A steady state solution corresponding to a zero-current state can be obtained 

A. 

dv 



by putting the term of ^ to 0. By solving eq. (H) we get a stationary 



distribution, 

p{v) oc {D^ + ny)^-^, (5) 
which has power law tails for large v, 

p{v) oc \v\''/^^-\ (6) 

For convenience of notation we define an exponent of the power law tail as 

p{v) oc \v\-^-\ (7) 

Comparing eq. (||) with eq. (0) yields 

P = -f. (8) 

It is useful to introduce a cumulative distribution function (CDF) defined as 

/ — \v\ roo 
p(v')dv'+ / p(v')dv'. (9) 



Then the CDF of v in eq. (|1|) is given as 

P(> \v\) oc \v\-^, (10) 

namely, the slope in log-log plots directly gives the power law exponent (3. 
The corresponding PDF of eq. is given as 

p(„)oc(l4)-^'^-'^ (11) 

where s = \f^- Namely p{v) has power law tails for v > s. 

3 The circuit 

From the assumption for the multiplicative noise in the above theory it is 
important that z/(t) takes both positive and negative values to realize the 
power law tails. In an electrical circuit this means that it is necessary for 
the circuit to include both positive and negative feedbacks. We solve this 
problem by using an analog multiplier. A block diagram of our circuit and an 
implemented noise generator are shown in fig.|I|. In the figures Vo represents 
the output voltage, and is the output directly from the noise generator. 
As seen from this figure the circuit contains the noise generator, an analog 
multiplier (Analog devices, lOMHz, 4-quadrant) and an operational amplifier 
(National Semiconductor, LF157) for integrator. 

The output of the noise generator plays a role of the multiplicative noise 
in eq. (|I]). In the noise generator a shot noise between the zener diode in fig. 
|I| is amplified by an operational amplifier, it passes through a high pass filter 
to be taken out from the circuit. LF157 has the 20M product of a voltage 
gain (G) and a bandwidth (B). The bandwidth is given by B=200kHz in 
the noise generator because we put G=100. Thus we expect a frequency 
characteristics of the noise generator to be up to 200kHz . We perform the 
Dositive and negative feedback bv multiulvine v„ with the outDut of the noise 



negative input with the operational amphfier for integrator. Although the 
additive noise term is not explicitly added in the circuit, it derives from 
either thermal noises of the operational amplifier or from an external electro- 
magnetic noises. 

An equation equivalent to the circuit diagram of fig. |l] is given as 

where is a factor of the multiplier with k = 1/10, and C,{t) represents the 
additive noise effect. The strength of multiplicative noise fi(t) depends on 
the value of a variable resistor in fig. |l] because is a factor of fi(t) in 
eq. (|l^). Therefore, we expect that the power law exponent for the output 
is a function of i?^ like eq. @. For the output of the noise generator /i(t) we 
show a PDF and an autocorrelation coefficient in figs. || and ^ respectively. 
The autocorrelation coefficient is defined as 

,,r , , (Mf + T)Mi))-(Mt + T))(A'W> ,,,, 

= wm - mr ' 

We sampled outputs of the noise generator through 12bit AD-converter (Mi- 
croscience, ADM-652AT) and processed them in a computer. A sampling 
frequency is 125kHz. We find that the PDF can be approximated by the 
Gaussian distribution and that the autocorrelation coefficient decays quickly. 
This result means that the frequency characteristics of the noise generator is 
up to around 200kHz. Therefore we consider as a Gaussian white noise 
under 200kHz. 

4 Results and Discussions 



We observe the output of the circuit Vo in respect of Ry since we expect 
that the power law exponent /? is a function of R^ from the theory. We 
detected the outDut throueh a 12-bit AD converter and Drocessed it as dieital 



we explain the method of adjusting an offset of the operational amplifier for 
integrator. We adjust the variable resistor added to the operational amplifier 
for a temporal average of output (v) to be almost likely zero when the value 
of i?„ is the smallest. As shown in fig. ^ the reason is that the temporal 
average shows exponential decay for Ry. We keep the zero-average of v{t) 
throughout all observations in this way. 

We show examples of temporal output at -R^, = 5 and = 150 in fig. 
^. We find that amplitude of the output v(t) depends on the value of Ry. 
The PDF for the temporal output on Ry = 5Q and Ry = 150Q are shown 
in fig.^ The tail's form of the PDF depends on the value of Ry. To see the 
difference between the positive and negative tails in detail we show log-log 
plots of the PDF in the positive and negative domains in fig. ^ It is obvious 
that the slope of the tail in the positive domain looks different from that in 
the negative domain. There are tow possible reasons on this effect. One is 
that an average of the additive noise is not completely zero. The other is 
that a distribution of the additive noise is asymmetrical. 

To roughly estimate the dependence of the exponent P on Ry we show 
log-log plots of CDF of v{t) at Ry = 5,25,50, 100, 150fi in fig. |. Log-log 
plots of CDF have a liner part for about a decade. We clearly find that the 
exponent /3 depends on Ry as expected qualitatively. However, we can not 
apply eq. (|) directly to these experimental results because the multiplicative 
noise generated by the noise generator is not an ideal white noise but is 
correlated at some frequencies. The authors have proved a relation to the 
power law exponent of RMP with correlated multiplicative noise in a discrete 
time version |]TT|], yet, in the continuous time version the theoretical relation 
between the exponent and the correlated multiplicative noise is an interesting 
open problem. 

Frnm thfi bpndinp" nnints nf CDFs in fiprlSl wp, ran fistima.tfi .s ke: 1 fnr 



eq. Namely, the strength of the additive noise is 10^^ times as large as 

that of the multiplicative noise. 

The autocorrelation coefficient of outputs, defined in eq. ( |T3| ) , are shown 
in fig. ^. A correlated time of the autocorrelation coefficient Ki{t) on 
Ry = 5Q is longer than the case of = 150f2. This implies that the 
correlation time depends on the value of R^. In economics it is often more 
interesting to observe the autocorrelation coefficient of squared outputs than 
the standard autocorrelation coefficient |jl2| . The autocorrelation coefficient 



of squared outputs called the volatility correlation is defined as 

, . {v{t + T)Mtf)~{v{t + Tf){v{tf} 

"^^^^ = (.(*)') - m^-r ■ 

As shown in fig. [1^ the correlation time of the autocorrelation coefficient 
K2{t) also depend on the power law exponent. The reason why Ki{t) or 
7^2 (t) have longer correlation time for smaller value of Ry can be simply 
answered in the following way. Comparing the time series in fig. |^ a and b we 
find that the motions are basically characterized by repeating the processes 
of starting from nearby zero, reaching the peak voltage and returning to 
zero again. As the maximum peak levels for R^ = 5Q are roughly several 
times larger than the case of Ry = 150f2, each process of going up and down 
roughly takes several times longer time leading the longer correlation time. 

5 Conclusion 

We proposed an electrical analog circuit generating the fluctuation of which 
a probability density function has the power law tails. The circuit includes 
a noise generator, an analog multiplier and an operational amplifler for in- 
tegrator. In our proposal circuit fluctuations with an arbitrary exponent of 
the power law can be obtained by tuning the value of resistance. We ob- 
served temDoral outDuts of the circuit and calculated the orobabihtv densitv 



the squared outputs. We show the probabihty density function with power 
law tails. The probability density function has different tails in positive and 
negative domains because the average of the additive noise is non-zero or its 
distribution function is asymmetrical. We conclude that the correlation time 
of the autocorrelation coefficient and that of the squared outputs depend on 
the value of resistance in the circuit. We expect that the proposed circuit is 
applicable to generate fluctuations having power law distribution in a much 
cheaper way than any digital computing methods. Moreover fluctuations of 
our circuit may be of use for risk estimation in foreign exchange or stock 
market in the near future. 
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Figure 1: Block diagram of all circuits with a noise generator. The circuit 
contains a noise generator, a analog multiplier and an operational amplifier 
for integrator. Rf — lOOkf), C = 10pF,-Rt, = 20017. A variable resistor under 
the operational amplifier is for adjustment of the offset. 



Figure 2: Semi-log plots of the probability density function of the multiplica- 
tive noise ii{t). Filled circles represent observation. A solid curve represents 
a Gaussian distribution with the same average and deviation as observa- 
tion. From numerical estimation the average is 0.005103, and the variance 
0.006846. 



Figure 3: Autocorrelation coefficient of the multiplicative noise //(t) for 
125kHz sampling frequency. Filled circles represent observation. 



Figure 4: Semi-log plots of relation between a temporal average of the output 
(v) and the value of Filled and Unfilled circles represent observation with 
different offset voltage in the operational amplifier. A solid line represents 
(v) = 4.5 exp(— 0.012i?^). The fit of the data to exponential is good from 
to 150Q. 



Figure 5: Signal outputs for = 5Q (a) and Ry — 150Q (b). The amplitude 
for Ry — 50, is larger than for i?„ = ISOQ. 



Figure 6: Semi-log plots of probability density functions R^ = 5Q (unfilled 
circles) and Ry = 150 fl (filled circles). The probability density function for 
Ry — 5Q has wider tails than for i?„ = 150Q. 



Figure 7: Log-log plots of the probabihty density function for R^ — 5Q. 
Unfilled circles represent it in the positive domain and filled circles in the 
negative. 



Figure 8: Log-log plots of CDF of the output Vo{t) at i?„ = 
5, 25, 50, 100, 150n. CDFs at Ry = 50, 100 show straight lines with differ- 
ent slopes between a decade. 



Figure 9: Log-log plots of autocorrelation coefficients for = 5Q (unfilled 
circle) and Ry = 150fl (filled circle). The autocorrelation coefficient for 
Ry — 5fl has a longer decay time than for Ry — 1500. 



Figure 10: Log-log plots of autocorrelation coefficients for squared output. 
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